New Class of Gravitational Wave Templates for Inspiralling Compact Binaries 
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Compact binaries inspiralling along quasi-circular orbits are the most plausible gravitational wave 
(GW) sources for the operational, planned and proposed laser interferometers. We provide new class 
of restricted post-Newtonian accurate GW templates for non-spinning compact binaries inspiralling 
along PN accurate quasi-circular orbits. Arguments based on data analysis, theoretical and astro- 
physical considerations are invoked to show why these time-domain Taylor approximants should be 
interesting to various GW data analysis communities. [? ] 
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Introduction. — Inspiralling comparable mass com- 
pact binaries are the most plausible sources of gravita- 
tional radiation for the operational, planned and pro- 
posed laser interferometric GW interferometers. GW 
data analysts, analyzing noisy data from the interferom- 
eters, require accurate and efficient temporally evolving 
GW polarizations, h + (t) and h x {t), the so-called GW 
search templates. It is expected that weak GW signals, 
buried in the noisy interferometric data, should be ex- 
tracted by employing the technique of 'matched filter- 
ing'. This is an optimal technique if and only if one 
can construct search templates that accurately model ex- 
pected GW signals from astrophysical sources, especially 
in their phase evolution. Till the late stages of binary 
inspiral, GW signals may be accurately modeled using 
the post-Newtonian (PN) approximation to general rel- 
ativity. The PN approximation to the dynamics of in- 
spiralling compact binaries, usually modeled to consist 
of point masses, provides, for example, the equations of 
motion as corrections to the Newtonian one in terms of 
(v/c) 2 ~ Gm/c 2 r, where v, to, and r are the character- 
istic orbital velocity, the total mass, and the typical or- 
bital separation, respectively. In PN computations, it is 
customary to treat a non-spinning inspiralling compact 
binary to consist of two point masses moving in quasi- 
circular orbits. These PN computations, to date provided 
four quantities that are required to do astrophysics with 
GW interferometers. For inspiralling compact binaries, 
the relevant four quantities are the 3PN accurate dy- 
namical (orbital) energy £(x), expressed as a PN series 

in terms of x = (Gwi^pn/c 3 ) 2 3 , o>3pn(£) being the 
3PN accurate orbital angular frequency, the 3.5PN ac- 
curate expression for GW energy luminosity C(x) and 
the 2.5PN amplitude corrected expressions for h+(t) and 
h x (t), written in terms of the orbital phase cj> and x 

GW data analysts employ these inputs to construct 



various types of search templates and let us take a closer 
look at the so-called TaylorTl and TaylorT2 waveforms 
implemented in the LSC Algorithms Library (LAL) 0]. 
These two template families employ the following expres- 
sion for the so-called restricted PN waveform 
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where the proportionality constant may be set to unity 
for non-spinning compact binaries. At a given PN order, 
the above mentioned two families provide two slightly 
different ways to compute w(t) and <fi(t). The TaylorTl 
family numerically solves the following two differential 
equations: 

dd>(t) , . du)(i) „, , / d£ 

J ¥ = "^ -^ = -«/^' (2) 

where, for example, C(oj) and £ are respectively the 
3.5PN accurate GW energy luminosity and the 3PN ac- 
curate orbital energy for TaylorTl 3.5PN waveforms. In 
other words, for a given PN member of the TaylorTl fam- 
ily, uj(t) and <f>(t) are computed by numerically solving the 
related approximants in Eq. ((2). To construct a mem- 
ber of TaylorT2 family, say TaylorT2 3.5PN, we require 
3.5PN (Taylor expanded) accurate version of dcu(t)/dt, 
appearing in Eq. j2J . The differential equations that de- 
fine u>(t) and (j)(t) for TaylorT2 3.5PN waveforms can be 
symbolically displayed as 



d<t>(t) du(t) 
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dt v ' dt 

+ 0(v) + (D(v 3/2 ) + 0{v 2 ) + 0{v b ' 2 ) 

+ 0(is 3 ) + 0(is 7/2 ) 



11 This is roughly the version submiited to Phys. Rev. Lett, on 
18/08/2007. A substantially modified version that also addresses 
objections of the referee is under preparation. 
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where v — 1/c 2 is a PN ordering parameter and the ex- 
plicit expressions for these PN contributions may be ex- 
tracted from Refs. [l|. In the above equation, the chirp 
mass A4 = to ry 3 ^ 5 , where r\ is the usual symmetric mass 
ratio and to being the total mass of the binary. 
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In this paper, we provide prescriptions to compute 
three new types of time-domain Taylor approximants that 
should be, in our opinion, interesting to various GW data 
analysis communities. Let us first list the salient fea- 
tures of these new templates that also employ an expres- 
sion similar to Eq. (fT]) to generate waveforms. The three 
important features of our Taylor approximants are the 
following. The first point is that, in comparison with 
TaylorTl and Taylor T2 waveforms, for a given GW fre- 
quency window and at a given PN order, our prescrip- 
tions will provide more accumulated GW cycles. Further, 
our approaches to compute h(t) are numerically as cheap 
(expensive) as TaylorTl and Taylor T2 waveforms. Let 
us consider the second point. It is desirable to construct 
GW templates using the mathematical formulation em- 
ployed to construct the (heavily employed) PN accurate 
relativistic D amour- Deruelle timing formula for binary 
pulsars [H . This is because formally GW phasing for in- 
spiralling compact binaries and timing of relativistic bi- 
nary pulsars are quite similar. Our construction of these 
new Taylor approximants are indeed influenced by the 
GW phasing formalism, available in Ref. Q, that pro- 
vided a method to construct GW templates for compact 
binaries of arbitrary mass ratio moving in inspiralling ec- 
centric orbits. We recall that the techniques adapted in 
Ref. Q were influenced by the mathematical formulation, 
developed in Ref. [5|], to compute the D amour- Deruelle 
timing formula. Finally, a recent preliminary investi- 
gation indicates that our new Taylor approximants, at 
the dominant radiation reaction order, should be very 
efficient in capturing GWs from compact binaries inspi- 
ralling along PN accurate and mildly eccentric orbits [6j . 
This is, in our opinion, a very attractive feature for GW 
data analysts as GWs from inspiralling (astrophysical) 
compact binaries should have some tiny eccentricities, 
when their GWs enter the bandwidth of laser interfer- 
ometers. 

Let us describe how we construct these new types 
of PN accurate time-domain Taylor approximant GW 
search templates. GWs from inspiralling (astrophysi- 
cal) compact binaries will have some tiny eccentricities 
around orbital frequencies of 20Hz. For example, using 
Ref. it is not that difficult to show that the orbital 
eccentricity of the Hulse- Taylor binary pulsar when its 
orbital frequency reaches around 20 Hz will be ~ 10~ 6 . 
Therefore, let us take a closer look at how one can de- 
scribe, in a PN accurate manner, eccentric orbits, mo- 
tivated by the fact that GW phasing requires accurate 
orbital description. Couple of decades ago, it was demon- 
strated that associated with a PN accurate non-circular 
orbit, there exists two gauge invariant quantities, if ex- 
pressed in terms the conserved orbital energy and an- 
gular momentum of the binary 0]. These are the PN 
accurate mean motion n and k that measures the ad- 
vance of periastron in the time interval T, T being the 
radial orbital period such that n — 2ir/T. It is quite 



convenient to define these quantities in the PN accu- 
rate Keplerian type parametric solution to the conserva- 
tive PN accurate compact binary dynamics, available in 
Refs. 8]. When eccentricity parameter, say time eccen- 
tricity e t , associated with the PN accurate Keplerian type 
parametrization approaches zero, one can define PN ac- 
curate orbital angular frequency oj = dcj>/dt = n(l + k). 
This implies that for PN accurate circular orbits, the 
angular part of the orbital motion is simply given by 
<f> — (bo = n x (1 + k) x (t — to). To 3PN order, using 
Ref. 0, in the limit et — > 0, we have 
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where £ = Gmn/c 3 . Let us now compute employing 
PN accurate expressions for £ (x) and C(x), available in 
Refs. [H, the following 3PN accurate expression for the 
orbital energy £ and 3.5PN accurate GW energy lumi- 
nosity C, in terms of £, as 
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where £ = —2E, E being the dimensionless non- 
relativistic energy per unit reduced mass [J] and 7 be- 
ing the Euler's gamma. We are now in a position to 
construct, in our terminology, TaylorKl 3.5PN and Tay- 
lorK2 3.5PN restricted PN waveforms. In our approach, 
the form of the restricted PN waveform, Eq. (fT]), becomes 

I Gm n(t ) \ 2 ^ 

h{t) oc ( "3 ' J cos 2 This is allowed because 

at Newtonian order lo = n and the amplitude is indeed 
Newtonian accurate in Eq. (fT]). For TaylorKl 3.5PN ac- 
curate waveform, nit) and (b{t) are numerically obtained 
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using the following two differential equations 



dt U V dn 



(6a) 
(6b) 



To construct our TaylorK2 3.5PN waveforms, as ex- 
pected, we Taylor expand, in terms of £, the RHS of 
Eq. (|6bj) and this leads to 
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Let us now specify, for example, the limits of integration 
for n to construct TaylorKl 3.5PN and TaylorK2 3.5PN 
waveforms. For initial LIGO, it is customary to use oj, 
and u)f, the initial and final final values of u>, to be 40 7r 
Hz and (6( 3 / 2 ) m) _1 Hz, where u>f is twice the conven- 
tional orbital angular frequency of the innermost stable 
circular orbit for a test particle around a Schwarzschild 
black hole. With these inputs, the initial and final values 
of n, denoted by n-i and n/, are numerically computed us- 
ing Eq. ffl. This is justified because of the observation 
in Ref. [lfj that the quadrupolar GW frequency from a 
compact binary, having PN accurate orbital motion, ap- 
pears at (1 + fc) n/n. At 3PN order, for a compact binary 
having m = 11.4Mq and 77 ~ 0.108 we have rii ~ 111.32 
Hz and 77./ ~ 679.3Hz. In our approaches to construct, 
for example, TaylorKl 2PN and TaylorK2 2PN wave- 
forms, we use only the 2PN accurate relation connecting 
to and n. 

Let us now compute in the time domain the accumu- 
lated number of GW cycles, Afcw, m a given GW fre- 
quency window, by numerically integrating Eqs. ^ and 
([7]) representing temporal evolutions for TaylorKl and 
TaylorK2 waveforms at four different PN orders, namely 
2PN, 2.5PN, 3PN and 3.5PN orders, for three canonical 
compact binaries usually considered in the GW litera- 
ture [we restrict these orbital evolutions such that emit- 
ted GWs are in the GW frequency window defined by 
40 Hz and (6 3 / 2 7rm) _1 Hz]. Let us also compare these 



Ngw with what is expected from TaylorTl and Tay- 
lorT2 waveforms at these four different PN orders . The 
numbers, relevant for initial LIGO, are listed in Table U 
where we compare New resulting from TaylorK2 and 
TaylorT2 prescriptions [results are similar while compar- 
ing TaylorKl and TaylorTl] 



TABLE I: Accumulated number of GW cycles, relevant for 
initial LIGO, for three types of canonical binaries at four dif- 
ferent PN orders using TaylorK2 and TaylorT2 waveforms. 
The values of Mgw arising from TaylorT2 waveforms are 
given in parentheses. We note that TaylorK2 waveforms pro- 
vide more Mgw compared to TaylorT2 waveforms. For high 
mass binaries, the convergence of Ngw is not that pronounced 
for TaylorK2 waveforms compared to TaylorT2 waveforms. 



mi/M© : m 2 /M Q 


1.4 : 1.4 


1.4 : 10 


10 : 10 


2PN 
2.5PN 
3PN 
3.5PN 


1616.4 (1613.5) 

1613.5 (1605.8) 
1623.4 (1616) 

1620.6 (1615.4) 


345.6 (333.8) 
333.8 (333.1) 
347.3 ( 330.9) 
342.4 (330.5) 


57 (52.6) 
53.8 ( 52.6) 
57.6 (52.9) 
56.2 (52.5) 



We are aware that LAL also provides routines to create 
TaylorT3 waveforms. In this prescription, both 4>(t) and 
ui(t), appearing in Eq. ([1]), are given as explicit PN accu- 
rate functions of time. These explicit time dependencies 
are usually expressed in terms of the so-called 'adimen- 
sional' time variable 9 = 5 c G v m (t c — t) , where t c is the PN 
accurate coalescence time. It is indeed possible for us to 
compute n(t), using Eqs. |J7J|, as a PN series in terms of 
6. However, we are reluctant to repeat what is done in 
TaylorT3 waveforms to get <p(t) with the help of Eqs. 0. 
Observe that radiation reaction and hence temporal evo- 
lution of 77 first appears at 2.5PN order and therefore, in 
our opinion, it is better to keep dcf>/dt to at least 2PN or- 
der in Eqs. (J7J) to be consistent in a PN way[see Refs. 
where similar approaches are employed]. 

It is important to note, while constructing these time- 
domain Taylor waveforms, that we employed the follow- 
ing two arguments. The first one is the standard argu- 
ment that equates the rate of decrease of the conserved 
orbital energy of a compact binary to the opposite of GW 
luminosity. However, for constructing TaylorTl, Tay- 
lorT2, TaylorKl and TaylorK2 waveforms, one requires 
additional PN accurate relations relating u> (or n as the 
case may be) to the conserved orbital energy. Further, 
we speculate that the two different ways of computing 
dw/dt, enforced in TaylorTl and TaylorT2 waveforms, 
may be based on the fact that observationally dui /dt (or 
the above mentioned standard argument) is only tested 
to the Newtonian radiation reaction order by the accurate 
timing of binary pulsars. Therefore, it is natural to ask if 
we can construct h(t) employing only the energy balance 
argument. This is indeed possible as demonstrated below 
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where t = tc 3 /Gm and C = £■ We call the resulting h(t) 
as TaylorEt waveforms. The values of £ corresponding 
to uji and LOj can numerically evaluated using the RHS 
of Eq. (gb| for d(j>/di = Co. 

We evaluated Now associated with TaylorEt 3.5PN 
waveforms for the three canonical compact binaries and 
the numbers are the following. For neutron star binaries, 
m = 2.8M Q and n = 0.25, M GW = 1617.4 and for the 
usual black hole-neutron star binaries, m = 11.4M and 
r] = 0.108, we have Mow = 335.4. For typical stellar 
mass black hole binaries, m — 20M Q and r\ = 0.25, one 
gets Ngw = 54.0. It is interesting to note that we get 
larger Maw, compared to TaylorT 3.5PN waveforms and 
lower values compared to TaylorK 3.5PN waveforms. It 
should be related to the fact that it takes more time 
for TaylorEt prescription to reach ujf form LOi compared 
to TaylorTl (or TaylorT2) approach and the opposite 
is true for the cases of TaylorKl (or TaylorK2). The 
observation that TaylorEt waveforms also provide more 
number of GW cycles in a give GW frequency window, in 
our opinion, makes it our third prescription to compute 
/'('). 

Conclusions. — We provided new ways of construct- 
ing restricted time-domain PN accurate waveforms for 
non-spinning compact binaries inspiralling along PN ac- 
curate quasi-circular orbits. Our prescriptions employed 
PN accurate expressions for the conserved orbital energy 
and GW luminosity, available in Refs. [l| , in a democratic 
manner and heavily depended on certain PN accurate 
gauge invariant quantities, first introduced in Ref. Q- 
These template waveforms provide more number of ac- 



cumulated GW cycles in a given GW frequency window 
and may be useful in detecting GWs from inspiraling 
compact binaries that should have 'teeny-weeny' orbital 
eccentricities. Further, our approaches are influenced by 
the way PN accurate D amour- Deruelle timing formula 
was constructed. Therefore, we feel that our TaylorKl, 
TaylorK2 and TaylorEt waveforms should be of certain 
interest to the practitioners of LAL. Further, we feel that 
our restricted PN waveforms should be useful for the the 
recently initiated mock LISA data challenge task force. 

The data analysis implications of these templates, rele- 
vant for both ground and space based GW detectors, are 
under active investigations in collaborations with Stas 
Babak, Sukanta Bose, Christian Rover and Manuel Tcss- 
mer. The GW phase evolution under our prescription is 
also being compared with its counterpart in numerical 
relativity based binary black inspiral. 
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